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Abstract 

In this paper, we theoretically investigate how laser fields change the nature of topological Kondo 
insulators(TKIs). By employing a prototypical model of TKIs, we treat the effect of the laser 
fields with Floquet theory, which gives effective description under high frequency laser fields. 
We derive the effective model of TKIs under the laser irradiation and discuss its topological 
properties. We demonstrate a possible realization of Floquet Chern insulators specified by 
various values of Chern number and reveal how the topological phase changes with increasing 
the laser light intensity. Furthermore, it is shown that Floquet Weyl semimetals, which have 
some pairs of Weyl nodes protected topologically, can emerge in the three-dimensional case. 
We explain how the Weyl nodes are created with varying the strength of the laser field. 

Keywords: photo-induced phase transition, topological matter, topological Kondo insulator, Weyl 
semimetal 


1 Introduction 

In recent research for new topological phases of matter, topological Kondo insulators (TKIs) 
have been proposed and have attracted lots of research interest 0. Unlike the other well studied 
topological insulators, e.g. Bi 2 Se 3 , TKIs are strongly correlated electron systems that have 
robust conducting surface states protected by time-reversal symmetry. Theoretical study by 
Dzero et al. [2j opened a novel direction of research for interacting topological matter. From the 
experimental viewpoint, SmB 6 is believed as a strong candidate for TKIs. There are many pieces 
of evidence supporting that SrnBg has topologically robust surface states at low temperatures, 
while negative results are also obtained experimentally. Although it is still an important open 
question to figure out whether the Kondo insulator SinBg is topological or not, TKIs certainly 
provide a versatile platform to study topological phases of matter with strong correlation. 

Besides intensive studies on topological properties in equilibrium, in recent years, there 
has been tremendous progress in emergent quantum phases of matter induced by laser light 
applications. Especially, it has been intensively investigated what topological phases are induced 
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by laser light. For example, it is known that the application of linearly polarized laser can induce 
a topological state in a semiconductor quantum well, which is initially in the trivial phase j7j. 
It exemplifies a realization of nonequilibrium states corresponding to Z 2 topological insulators, 
called Floquet topological insulators. There are also other examples showing that the system 
can be driven to nonequilibrium topological phases in graphene [6] or cold atomic systems [8]. 

In this study, we theoretically investigate the topological nature of TKIs under laser irra¬ 
diation. TKIs have the unique hybridization term which depends on the (pesudo)spin and the 
momentum of electrons. It originates from the spin-orbit interaction(SOI). Then a new kind 
of phenomenon can be induced by the cooperation of the effect of laser light and SOI in the 
laser-irradiated TKIs. The paper is organized as follows. After introducing the model of TKIs 
and the method to treat the effect of high frequency laser irradiation, we derive the effective 
model and discuss its topological properties. We address a realization of the laser-induced 
topological states of matter such as Chern insulators and Weyl semimetals, which are referred 
to as Floquet Chern insulators and Floquet Weyl semimetals |l'3j. 


2 Model and Method 


2.1 Model Hamiltonian of topological Kondo insulators 


In order to discuss the laser-induced effects on TKIs , we introduce a simple model{TD], which 
consists of the conduction d-electron band and the localized /-electron band. This model is a 
variant of the periodic Anderson model of which hybridization terms have nontrivial depen¬ 
dences on the (pseudo)spin and the momentum due to SOI in /-orbit. Assuming the cubic 
lattice structure similar to SmBg, which is a strong candidate material for TKIs, the model 
Hamiltonian reads 


%tki = 


fee b.z. 


£d(fe) 1 
V(°)*(fc) • 


y(°)(fe) ■ c 



'it 


(i) 


with 


ed(k) = — 2 td(cosk x + cos k y + cos k z ), (2) 

e^\k) = — 2tf(cosk x + cos k y + cos k z ), (3) 

= iV^°\ai sinfc x , ai sink y , 0,3 sin k z ), (4) 


where 1 and <x are 2x2 identity matrix and Pauli matrices. We use (d*., f *.) as a short hand 
notation for (dfc-j-, dfcj.,/fc-j-,/fcj.), which are annihilation operators in the conduction and the 
localized band respectively. Here / and / stand for the (pseudo)spin. The fact that e<j(fc) and 
ef(k) are parity even and V ( k ) is parity odd guarantees time reversal symmetry of this system 
and we can define a Z 2 topological index. 

We incorporate the correlation effect in /-orbit with a Gutzwiller-type mean field treatment, 
following Dzero et al. [2|. We thus assume a renormalized band structure, which results in the 
mean-field Hamiltonian, 


Hmf = E &(k)H M F(k)*(k) 

ke B.z. 

V f d i\ T ( ed(fe) 1 V(k)-a 
[ft) \V*(k)-a e f (k) 1 



( 5 ) 

( 6 ) 
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where 4'(ft) = (<4 t , d ki , / fcT , V = FK (0) = iV(a 1 smk x ,a 2 smk y ,a 3 smk z ) and e f (k) = 
( k ) + A. Note that & is a Gutzwiller-type renormalization factor and A represents an energy 
shift of /-orbital level by the interaction effect. 

We set the parameters of the model as tf = —0.25 td in this study and vary the parameters 
V, ef and a = ( 01 , 02 ,( 13 ) according to the situations. Especially we consider two cases of the 
hybridization parameter a in this study: a, 2 D = (2,2,0), a 3 D = (2,2, —4). Here a. 2 D(with fixing 
k z = 0) and a 3 o correspond to a two-dimensional (2D) TKI and a three-dimensional (3D) TKI 
respectively. 


2.2 Effective Hamiltonian in Floquet theory 

It is known that Floquet theorem is useful for analyzing the system which has time periodic 
Hamiltonian[5]. Floquet theorem is, so to speak, “Bloch’s theorem for time direction”: If the 
Hamiltonian is time periodic H(i) = H(f + T), the eigenfunction can be written by a product 
of an exponential function e _let and a time periodic function u(t). e is called “pesudo energy” 
that is defined in the range of —7 t/T < e < 1 t/T = w/2. Then, we can define the effective 
Hamiltonian in Floquet theory, of which eigenvalues are pesudo energy, as 


n oS = l -\ogU(T), (7) 

where the time-evolution operator U{t) is written as Texp * f* ?{(s)dsj . By definition, the 
effective Hamiltonian has only the information of t = 0, T, 2T, • • • , nT, ■ ■ ■ and gives “strobo¬ 
scopic description” of this system. In this description, the information of the time range of 
nT < t < (n+ 1 )T is neglected. However, if the time period is short enough, this Hamiltonian 
tells us the effective behavior of this system; namely it describes the asymptotic behavior in the 
high frequency limit. Then we can understand the nonequilibrium states in the high frequency 
limit with the effective equilibrium model. 

There are some methods to calculate the effective Hamiltonian. We adopt a perturbative 
approach following Kitagawa et aZ.pJj. With this approach, we can write down the effective 
Hamiltonian as 


'Hoff 


H 0 + 


OJ 


+ 0 



( 8 ) 


where H n = ^ / T ^y 2 dt'H(t)e mult . The second term in the case of the laser-irradiated systems 
can be understood as the second order perturbation process of the virtual photon absorption 
and emission in off-resonant light. We derive and investigate the effective model under the high 
frequency laser fields with this approach in the following section. 


3 Results and Discussion 

3.1 Derivation of the effective model 

We derive the effective model of TKIs under the laser field A = ( A x sinwf, A y sin(o;f + (f>), 0), 
where tf> is the polarization angle of the laser light. We focus on (f> = 7r/2(circularly polarized) 
case in this study. In this case, the circularly polarized laser light breaks the time-reversal 
symmetry of the original system. The symmetry class is changed by the laser irradiation 
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and thus it is expected that the application of the laser light makes the nature of the system 
dramatically different. The effect of the laser light is taken into account by minimal coupling 
k —> k — A, which corresponds to introducing a Peierls phase t l: j —>• t-ij in the real space 
representation. Using H(t) = Jfke bz ^(fe)’HMF(fe — (k) , we calculate the effective 

Hamiltonian with the method in Sec |2.2[ The explicit form of the effective Hamiltonian obtained 
for <f> = n /2 is 

„ _ V (4\ T ( (V(fc)+ *(*))• „ \ 

JiU/ ((v(*=) + *'(*))•"■ i/(*)l + »#K ) W. 


where 

=-2t d (J 0 {A x ) cos k x + J 0 (A y ) cos k y +cos k z ), (10) 

e/(fc) = ef — 2tf(Jo(A x ) cos fca; + Jo(A y ) cos k y + cos/c z ), (11) 

V(fe) = jU(aiJo(^)sin^,a 2 Jo04z) s in&y,a3sinfcz), (12) 

^ /n 2U 2 J\{A X )Ji(Ay) * *\ i i 

®B(k) = -—(aia 2 + 020 !) cosk x cosk y , (13) 

3>(fe) = ( ai cos fcj; sin —a 2 sin cos k yi 0), (14) 

UJ 

and J n (x ) is the n-th Bessel function. 


We can see some effects of the laser light in eq.Q. i) Photo-renormalization effect: We 
can see the renormalization of the hopping integrals td,tf an d the hybridization term V by 
the 0-th Bessel function. Regarding hopping integrals, this effect is nothing but “dynamical 
localization” JTj , which claims that the hopping integral of electrons is reduced by the laser light 
because the frequency is too high for the motion of electrons to follow the AC electric fields of 
the laser light. The effect on the hybridization term can also be understood as an analog of 
dynamical localization since TKIs have only inter-site hybridization unlike the Kondo insulators 
described by the conventional periodic Anderson model, ii) Photo-induced hybridization: From 
the commutator term in eq.([8]), we obtain new hybridization terms 3>(fc). They are derived 
from spin-dependent hybridization terms originating from SOI. The forms of these terms reflect 
the 7r/4-rotational symmetry ( k x —¥ k y , k y —> —k x ) since the laser light has the same symmetry, 
iii) Photo-induced “magnetic” field: 4'^(fc) are derived from the commutator of the x and y 
components of hybridization terms and play a role of the Zeeman coupling to (pesudo)spins. 
Then we can regard them as momentum-dependent “magnetic” fields originating from SOI. 
These terms break time-reversal symmetry in a similar way to usual magnetic fields. In the 
following sections, we will see that the effect iii) particularly plays a vital role in inducing 
topological phase transitions. 

3.2 Topological invariant 

We investigate the topological nature of the effective model derived in the previous section. 
In general, the topological nature is characterized by a topological invariant, which is robust 
against any local perturbations without breaking the relevant symmetry of the system, e.g. 
time reversal symmetry or inversion symmetry. For this purpose, we use the Chern number C, 
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which is a topological invariant defined in two-dimensional systems. The Chern number can be 
calculated as m 

C = f del d« y Tr[ e ^(G9 M G- 1 )(G9 1 ,G- 1 )(Ga p G- 1 )], (15) 

24n- J_ 00 J B Z 

where the Green function G _1 (e, k) = ie — H e s (fe) and /z, ix, p run through e, k x , k y . As for 2D 
systems, it is known that the Chern number can characterize integer quantum Hall states or 
Chern insulating states. Moreover, topological characterization by the Chern number can be 
applied to 3D systems. Weyl semimetals, which have an even number of topologically protected 
gapless nodes (Weyl nodes), can be also characterized by the Chern number, which is calculated 
on a 2D plane in the 3D Brillouin zone. If we choose two planes which hold a Weyl node between 
them, we obtain the different Chern numbers from each plane. This difference corresponds to 
the topological number which the Weyl node has. 

3.3 2D case : Floquet Chern insulator 

First, we study the 2D case in which we use the parameters a 2 D = (2, 2, 0) and set k z = 0. 
We can calculate the Chern number of this system since the system is 2D and gapped. For 
several values of |A| and e/, we numerically evaluate the non-trivial Chern numbers. From this 
calculation, we obtain the topological phase diagram of this system, which is shown in FigJT] 
This phase diagram shows that the system can be Floquet Chern insulator (Quantum Hall 
insulator) phases, which have non-trivial Chern numbers. It is seen that they have some chiral 
edge modes corresponding to the Chern numbers of each phase. There are intriguing phases 
that have high Chern numbers G > 1 and |G| chiral edge modes. 
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Figure 1: Topological phase diagram of our effective model in the 2D case. The parameter is 
taken as V = 2.0 td- We can realize various phases that have different Chern numbers with 
varying the laser intensity A or the energy level of /-orbit ef. 


Furthermore, we directly confirm the existence of the chiral edge modes by numerical calcu¬ 
lation of the band structure. We calculate the band structure in the open boundary condition 
(OBC) in ^-direction and the periodic boundary condition (PBC) in y-direction. By this cal¬ 
culation, we obtain the band structure of Floquet Chern insulator that has G = 4 in Fig(2]A.. 
We can see from the figure that there are eight modes on the edges in ^-direction (Four nearly 
degenerate edge modes are at k = 0 and the other four edge modes are at k = ±7r). Further 
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investigation of the eigenfunctions of the edge modes elucidates that the four edge modes which 
propagate in the same direction are localized in the same side of the edges. We thus confirm 
that there are four “chiral” edge modes on each side, as shown in Figj2j3. We can control the 
number of these edge modes by tuning the laser light intensity |A|. 



Figure 2: (A) Band structure of Floquet Chern insulator calculated in OBC. The parameters 
are taken as |A| = 0.3 ,e/ = —4.0 1^, and V = 3.01 ^. (B) Schematic picture which describes the 
propagation of the four chrial edge modes when the circularly polarized laser irradiates the 2D 
system. We assume that the system size is finite in ^-direction and infinite in y-direction. 


Here, we remark difference between our study and the previous researches. It is known 
that photo-induced quantum Hall effect happens in graphene under the circularly polarized 
laser light [5]. This result is similar to our work in that the circularly polarized laser enables 
to acquire the non-trivial Chern number and induces the chiral edge modes. However, in 
the previous work, the effect of breaking the time reversal symmetry is induced by photo- 
assisted hopping between the next-nearest neighbors, whereas the photo-induced “magnetic 
held” originating from SOI breaks the time reversal symmetry in our model. There is not SOI 
in graphene and thus our proposal is different from theirs in the way to break the time reversal 
symmetry. 

3.4 3D case : Floquet Weyl semimetal 

Next, we study the 3D case in which we use the parameters <230 = (2, 2, —4). We numerically 
calculate the bulk band structure using the effective Hamiltonian, and show how the band 
structure changes with increasing the laser intensity in FigJ^fcj, and k y are fixed to zero). We 
show the initial state before the laser irradiation in Fig[3jA, and the band structure when the 
laser intensity is finite in Fig{3^3 and|3p. Weyl nodes are first created at k = (0, 0,0) and next 
at k = (0,0, 7 r), and some pairs of Weyl nodes can be seen in Fig(3j3 and[3]C. These results show 
the possibility that a Floquet Weyl semimetal[T3] is realized in this system. However, there is 
another possibility that these nodes are made by accidental band touching. To figure out this 
point, we examine the Chern number for several values of k z and find that the Chern number 
indeed changes if the k x -k y plane (fixing k z ) goes over the Weyl nodes when k z varies. This 
confirms that these nodes are topologically protected, thus ensuring the system to be a Weyl 
semimetal. 

Weyl semimetals receive a lot of attention recently as an exotic gapless topological phase of 
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Figure 3: Band structures for k z in the 3D case (k x and k y are fixed to zero). They are calculated 
in PBC for some values of laser intensity when tf and V are respectively fixed to —6.0^ and 
2.0fd- (A),(B), and (C) correspond to |A| = 0,O.ltd, 0.2td respectively. (A) represents the band 
structure of a TKI. In (B) and (C), we can see the Weyl nodes, which suggest that a Floqeut 
Weyl semimetal is realized. 


matter. There are interesting phenomena theoretically predicted in Weyl semimetals, such as 
the anomalous Hall effectfTT] and the chiral magnetic effect[l]. Unfortunately, there are only few 
experimental observations [5] of Weyl semimetals, and therefore it is desired to experimentally 
realize Weyl semimetals. An interesting feature of our proposal to realize Weyl semimetals in a 
nonequilibrium state is that we can control the number or position of Weyl nodes with varying 
the laser light intensity (see in Fig|3jU). This may allow us to control and observe characteristic 
surface states inherent in Weyl semimetals, called Fermi arc [12] which connects the different 
Weyl nodes. 

4 Summary 

We have examined how laser fields change the nature of TKIs. By employing a prototypical 
model of TKIs, we have treated the effect of the laser field with the effective Hamiltonian 
in Floquet theory. We have derived the effective model of TKIs under the laser irradiation 
and have revealed the effects of laser light. There are the effects originating from the spin- 
and momentum-dependent hybridization terms that are unique to TKIs. Then we have shown 
that the effects of laser light induce Floquet Chern insulators with various values of the Chern 
number and have elucidated how the topological phase transitions occur. Moreover, we have 
addressed a possible realization of Floquet Weyl semimetals, which have some pairs of topolog¬ 
ically protected Weyl nodes. 

In this paper, we have not discussed the important effects due to the laser-irradiation, i.e. 
the effect of heating, which should be crucial to realize our proposal in real experiments. This 
issue is now under consideration. 

This work is supported by KAKENHI (Grants No. 25400366 , No. 14J01328, and No. 
15H05855). M.N. thanks JSPS for the support from a Research Fellowship for Young Scientists. 
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